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Abstract 


This report describes the calibration of the school radio telescope. Section 2 outlines the 
theory of radio telescope operation. By executing a Sun observation the unknown system 
temperature can be determined. With this parameter, the performance of the telescope can 


be estimated (section 3). The calibration data and initial observation results are given in the 
Appendix. 
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1 Introduction 


Why we need calibration? A very demonstrative motivation of calibration was found in 1: 


As you likely know, every telescope is unique. One result of the uniqueness of 
individual telescopes is the difficulty of directly comparing measurements from 
one telescope with those from another. That is, Telescope A may record 280 
counts for the peak of a given spectral line, while Telescope B may record only 
100 counts. This is further complicated by the fact that even measurements 
taken on a given telescope, at a given frequency, can change over time. These 
changes can be the result of changes in, e.g. the telescope system temperature, 
the telescope response, and/or the atmospheric conditions. This means that 
even if you only observe your object with Telescope A, but measure it repeatedly 
over a year, you may discover that your object’s peak count varies from 280 one 
day to 250 or 300 on another day, and so forth. You then need to understand 
whether the source emission is truly varying with time or if the differences are 
due to changes within the telescope and equipment. 


In order to compare measurements between two telescopes, or even between 
one telescope taken at different times, we need a universal measurement 
system. That is, we need to be able to state that 280 counts from Telescope A is 
equivalent to X counts from Telescope B or equivalent to Y counts from 
Telescope A at a different epoch. This is the process of data (or telescope) 
calibration, and the rest of this chapter will be devoted to presenting the various 
methods available, both observationally and theoretically, to calibrate data. 


The task for calibration is to determine how the measured voltage levels are related to the 
flux of the observed sources. 


2 Calculations for Radio Telescopes 


2.1 Relation between measured voltage level and flux 


For cost reasons the telescope shall be calibrated without a local noise reference source. 
The receiver itself is calibrated and displays the incoming signal voltage level U in dBuV 
units. The voltage level is referenced to a resistance R of 75 Ohms. Hence the linear power 
of the received signal can be calculated: 


U[dBuV } d U? 
1 _19 2 and P[IW]=— 
(1) Ufuvy}=10 » VE 


If the telescope is pointed to a radiation source the power P is measured. Pointed away from 
any source a base power level Py is measured mainly caused by the noise of the telescope 
parts. The ratio between the power on source and off source is referred to as Y-factor. It can 
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be computed from the on source (UdB) and off source (U0dB) signal level values expressed 
in dBuV as follows: 





UdB-U 0dB 
(2) Y=—=10 © 


0 


The Y-factor is related to the antenna temperature of the observed source Tant, the 
temperature of the cold sky T.oa and the system temperature T.ys in the following way: 


Y i + a 
Q256 
fe + Teil 
For our telescope operated at ambient temperature we can assume that Tsys is much larger 


than Too. Therefore: 


y Ls + de 
4 = ee 
(4) Y= 


Sys 
Tantcan be determined for a known flux S as follows: 


S: A, 
eff 
(5) sie = K 2 
where Ager is the effective antenna aperture or the effective area of the antenna. Due to 
imperfections of the shape of the antenna, the effective aperture is smaller than the 
geometric size of the antenna. K is the Boltzmann constant equal to 1.38e-23 J/K. 





From equation 2, 4 and 5 we get the relation between the measured power value and the 
flux of the source: 


T gee. 





Solved for S: 


K.-T 
oy s-(2-1) Ae 
Py Aey 





Equation 6 shows, that the sensitivity of the telescope depends on Ae# and Tsys. The 
sensitivity increases if a larger antenna is used or if the system noise temperature is 
reduced. 


To estimate the sensitivity of our telescope, it is necessary to determine Ac and Tsys. Since 
the antenna is relatively small and solid we expect no significant changes of Ae# over the 
time. T,y; depends on the ambient temperature. Therefore it must be measured before the 
observations for calibration purposes. 


2.2 Determinination of Aeff 


The effective antenna aperture Aeff can be calculated with the following equation, taken 
from 2: 


471A 

RV 

G is the isotropic gain of the antenna, given in some datasheets of the antenna. / is the 
wavelength, where the dish is operated. 


G-#? 
3) Ay = 
(9) Ay =" 


(8) Ge eff 








2.3  Determinination of Tsys 


With the knowledge of Aes and the observation of a radio source with known flux, the system 
temperature can be calculated. Some well known point sources like Cassiopeia are used as 
calibration sources, since there flux is very stable. Due to the limited sensitivity we want us 
the Sun as calibration source. In contrast to Ac#, Tsys will change over the time. Therefore the 
calculated value is valid for a short period only. Nevertheless it may be an interesting task for 
students to perform a calibration. 


Equation 6 solved for Ts; provides: 
S+ Ag 








(10) Le 


2.4 Discussion about resolution and integration time 


The resolution 11 or minimum measurable temperature increase 12 of a radio telescope can 
be calculated by the "Dicke" expression 3. 





(11) ae 
P NB-t 

qe 
(12) je 


min JB-t 


Where B is the bandwidth and 7 is the measurement integration time. The expressions 11 
and 12 show, that the sensitivity can be improved by increasing bandwidth and integration 
time. 


Because this school radio telescope is based on the AMA 301 receiver, there is now option 
to change the bandwidth and integration time. The receiver bandwidth is about 20MHz in the 
used Sat-TV mode. This is a high bandwidth compared with ham radio receivers used in 
other radio astronomy projects. In the datasheets of the receiver there is unfortunately no 
information about integration time and measurement of the signal strength. After warm-up 
the absolute accuracy of the measurement is specified with +/- 1 dBuV, the values are 
provided with 0.1 dBuV resolution. 


The accuracy of the power measurement can be increased by computing the mean value of 
a number of single measurements. It requires an accurate tracking of the object position on 
the sky, which can be realized by using the tracking feature. To observe a weak radio 
source, a high absolute pointing accuracy of the telescope is additionally needed. It must be 
ensured, that the telescope is pointed to the source, before the long term integration shows 
a small increasing of recorded signal level. This accuracy is not provided by any rotator unit 
in a reasonable prize range. 


3 Calibration and Evaluation of the Telescope 


3.1 Computation of Aeff 


The effective antenna aperture Aeff is calculated as described in section 2.2. Unfortunately 
the manufacturer of our dish does not provide a value for the isotropic antenna gain G. The 
data of a similar antenna by Hirschmann are used instead: 

G is 41.5 dBi at 10.95 GHz, which corresponds to a wavelength / of 2.7378e-2m. 


4 EB _ 108" (2.747) 
oe 4n 4n 
The geometric size of a 1.20m antenna is about 1.13m?. As expected, the effective antenna 


size is smaller. Therefore the efficiency factor of the antenna is 0.745. This value is in the 
expected range. 


= 0,842546m’ 





3.2 Calibration of Tsys with Sun observation 


The results of a Sun observation (see Appendix 5.2 and 5.5) and reference values for the 
radio flux of the Sun at several frequencies (see Appendix 5.1) are used with equation 10 to 
determine the current value of Toys. 


The observation of Sun gives the for example the following values: 


402.1 sfu 1sfu = 10000J = 10°22 W/m2/Hz 
See Appendix 5.1 and 5.5 
[See3.t 0 


0.8425m? See 3.1 


1.389 J/K Boltzmann constant 
'P  st—i<‘t«isésSC«id:*«CY 3.5 dV See Appendix 5.5 (2005-08-05 — 16:20 
[Po ———«&X| 44.0 dBuV See Appendix 5.5 (2005-08-05 — 16:20 





With this data we compute a system temperature T.,, of 310.0 K. 


3.3 Estimation of telescope sensitivity 


The increasing of voltage level measured at the receiver while pointing from the cold sky to a 
radio source with the flux S can be computed from equations 2 and 6: 





SYS 


Ay 
(13) Ub -UddB = 10a ; 


With the values for Tsy; and Aes computed before and the flux value of Cassiopeia A at 10 
GHz of 1000 Jansky, we can predict the signal level at the receiver. Pointing to Cassiopeia 
the signal level will be increased about 0.009 dBuV compared to the cold sky. This weak 
signal cannot be detected in a direct way, because it is below the resolution of the device of 
0.1dBuV. Averaging between numbers of single measures can increase the resolution, but 
this method is difficult to use due to the restricted pointing accuracy of the telescope (also 
discussed in section 2.4. The figure shows the estimated signal level at the receiver versus 
the flux value of the observation object. The graph depends on the system temperature Ts 
which varies on environmental properties such as temperature. 


Telescope sensitivity for Aeff =0.84255m2 and Tsys =300K 
12 Titian) © fhe. fo tt fF Prime) © t beri ae 
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Level increasing [dBuV] 
oD 
T 

















4 4 
-Moon 
2 2] 
— geil i i pt —— oo 
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Flux [Jansky] 


3.4 Validation of the results at Moon observation 


Using equation 7 with the observation data of the moon (see appendix 5.5) we compute the 
flux of the moon to about 1.0e+004 Jansky. The intensity of thermal radiation of the Moon 
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depends on the Moon phase 4. Accurate radio flux values of the Moon could not be found on 
the Internet and in publications. A raw value can be taken from a diagram published in 2. 
The published value of about 5.00e+4 Jansky is above the value observed with our 
telescope. 


Source of the difference between published and observed flux value can be: 
e The dish is not pointed exactly on the radiation source while measure the on source 
level (the maximum was not captured) 
e Terrestrial radio sources as telecommunication services, near buildings can influence 
the results at telescope location in a city area. 
e The diagram 2 does not provide an accurate reference value for the flux of the Moon. 
e Heating of the LNB while pointing to the Sun at days with clear weather 
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5 Appendix 


5.1. Reference data for Solar fluxes 


The Learmonth Observatory (Australia) provides daily solar flux data on the Internet: 
http://www.ips.gov.au/Solar/3/4/2 
The fluxes are measured at several frequencies from 245 MHz up to 15.4 GHz. Between the 
observed values some other frequencies are interpolated. The values represent the quite 
Sun. If bursts are detected this frequencies are not useful for interpolation. The flux is given 
in solar flux units (sfu). 
1sfu = 10000 Jansky = 10°? W/m?/Hz 


The data from Learmonth observatory are provided in the format shown below. 


Quiet Solar (IFLUX) 


Last updated 05 Aug 2005 07:30 UT 


IFLUX : Background Solar Radio Flux 


Station Date Time Status Freq OS i rmlausc Quality 
Learmonth 05/08/5 03:48 final 245 17 ? burst 
410 24 good 

610 44 good 

1415 63 good 

ZAG) 114 good 

4995 149 good 

8800 ZS good 

15400 ple, good 








Interpolated value for 1300MHz: 60. 
Interpolated value for 1540MHz: 68. 
Interpolated value for 1707MHz: 74. 
Interpolated value for 2300MHz: 98. 
Interpolated value for 2401MHz: 102. 
Interpolated value for 2790MHz: 115. 
Interpolated value for 5625MHz: 166. 
Interpolated value for 6000MHz: 176. 
Interpolated value for 8000MHz: 231. 
Interpolated value for 8200MHz: 236. 
Interpolated value for 10400MHz: 312.3 


CHRO TES eo: 


onrUOATUY 


For our frequency of 12.6GHz we estimate a flux of about 402.1sfu by linear interpolation of 
the data in Excel. 
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An alternative Internet resource for radio flux values of the Sun is: 


http://www.sec.noaa.gov/ftpdir/lists/radio/7day_rad.txt 


5.2 Sun observation data (sample) 





log_2003 03 15 _07.txt 


61 
60 
59 
58 
57 
56 
55 
54 
53 
52 
51 
50 
16:13:26 16:16:19 16:19:12 16:22:05 16:24:58 16:27:50 


Time [hh:mm:ss] 








ay 





S 

















Level [dBuV] 



































Sun Measurement, Freq.= 12600 MHz, SAT-IF= 2000 MHz, High Band, Horizontal Polar. 
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5.3. Moon observation data (sample) 























Level [dByV] 
& 
KR 
aN 


























43,8 























43,6 
10:30:43 10:33:36 10:36:29 10:39:22 10:42:14 10:45:07 10:48:00 10:50:53 10:53:46 
Time [hh:mm:ss] 














5.4 Matlab script for the computation of the parameters 


The m-File "calib3.m" contains: 

e the computation of Ts,; based on a Sun observation 

e a plot of the telescope sensitivity 

e estimation of Moon and Cas A signal values 
The script can be used with Matlab (commercial software for technical computing) or Octave, 
a free Matlab clone available under Gnu Public License (GPL), distributed with Linux. 


SRadioastronomie at Schools 
SComputation of System temperature 


SParameters of the telescope 





Aeff=0.842546; seffective antenna aperture m2 
K=1.38e-23; SBoltzmann J/K 
R=75; Sresistance (not used) 


SObservation results 


UdB=60.3; SLevel measured for Sun [dBuV] 

U0dB=50.8; SLevel measured for cold sky [dBuV] 
Ssun=400.e-22; Sreference for Solar flux from Internet 
UmdB=44.9; smeasured value for Moon dBuV 

UOmdB=43.9; smeasured value for cold sky near Moon dBuV 


% CALIBRATION WITH SUN 
$----Computing Y faxtor linear------ 


Y=10% ( (UdB-UOdB) /10) ; 
disp(['Ysun (lin) =' 


, 


num2str(Y)J]); 


° 


3----Computing System temperature 
tl=Y¥=1; 

t2=Ssun*Aeff/K; 

Tsys=t2/tl1; 


disp(['Tsys =', num2str(Tsys) ]); 


= 


G ml 


G 


ELI 





THI 

















=] 








ole 


ESTIMAT! SENSITIVITY OF TI! 
Scompute a sweep of flux values 
S_ja=1:10000:5000000; 


S=S_ja*le-26; 





Sestimated signal level in dBuV vs. 
dU_log=10*1log10 ((S*Aeff/ (K*Tsys) ) +1 


Splot the results 
plot(S_ja,dU_log); 
title(['1 
num2str(Tsys), 'K']); 
xlabel('Flux [Jansky]'); 
ylabel('Level increasing 
grid on; 





[dBuV]'); 


ole 


ESTIMATION OF FLUX VALU! 








E SCOP! 


[Telescope sensitivity for Aeff =', 


ES FOR CASSIOP!I 


flux 
\e 


num2str(Aeff), 'm2 and Tsys =', 





EIA 


Swhich voltage level gives Cassiopeia at the receiver? 


Scas=1000e-26; Stabulate 


d flux for cassiopeia 


dP2log_cas=10*1log10((Scas*Aeff/ (K*Tsys) ) +1); 





10GHz (dBuV) , num2str(dP2log_cas)]); 





disp(['Estimated level for CAS A at 
SVALIDATION WITH DATA FROM MOON OBS! 
Ym=10% ( (UmdB-UOmdB) /10) ; SY-factor 
t3 = K*Tsys/Aeff; 

t4 = Ym-1; 

Smoon = t4*t3/10e-26; 


disp(['Calculated Moon flux (Jansky 


ERVATION 
linear 


’ 


) 


num2str(Smoon) ]); 
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5.5 Observation Data 


Frequency Solar Flux | Computed | Computed Weather Comments 
[sfu] Flux of 
Moon 
[Jansky] 
2005-08-04 : : : 407.0 : Tsys to high. 


Maximum probably 
not captured 
14:58 | 12600(H Moon _}444 {43.7 ff 11099.6 i oe 
15:09 | 12600(H Pot 1844.5 Pp 
12600 = 3140 | St 


2005-08-05 1S: 15 | 12600(H 54, 0 13. 8 402, 1 259.2 Sunny, Power level jumped 
to low? 20 C during the 
measurement, value 
may be corrupted 


15:32_| 12600(H ‘Moon _}44.1 {49.2 ff |9776.0 Se 
16:11_| 12600(H | CT CC 8888.0 ness) 
as 
a 








/16:20 |12600(H) [Sun [535 [440 [4021 [3100 | 
[16:34 [42600(H) [Moon [445 [439 | | «(7528.9 
2005-08-08 | 15:00 |12600(H) [Sun 1533 1436 | 399.0 [2924 | | Sunnyand =| 
| 16:20 |12600(H) [Sun [53.7 | 44.0 [3990 [2924 | Clouds, 16C | 
| 16:33 | 12600(H) | Moon | 44.7 0) 44.000 8377.0 (= = = 
2005-08-09 390.7 {3527 | | Sunny 18 C_| Receiver Warm-up? _| 
(3176 | Clouds 18 C | Fulltransit_ 
12600(H) | Sun peo27 — | Clouds 18:0. «|| Omettmesstiement 
14: 8s | (3346 | 
14: as (3436 | 
14:50 _| i a 
| 2005-08-11 | 1 12600(H) | Sun (1420 (515 3907 =~] 3013 | Clouds 20C | 


Solar Flux values are used from the Learmonth observatory see section 5.1. 
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